A steady state network model with a 1/k scale-free degree distribution 
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Using a steady state process of node duplication and deletion we produce networks with 1/k scale- 
free degree distributions in the limit of vanishing connectance. This occurs even though there is no 
growth involved and inherent preferential attachment is counterbalanced by preferential detachment. 
The mean field evolution is considered and the 1/k law is verified under certain approximations. An 
ansatz for the degree distribution is proposed on the basis of symmetry considerations and is shown 
to coincide well with the simulation data. Distributional forms other than power law are also shown 
to arise when the condition of perfect duplication is relaxed. 

PACS numbers: 89.75Fb, 89.75Hc, 87.23Kg 
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The frequent appearance of power-law degree distri- 
butions in real system networks is an intriguing phe- 
nomenon. It has been suggested that the scale-free na- 
ture of such inter-connected systems offers a level of ro- 
bustness and stability to the dynamics P,0 so is favoured 
in a selective sense. Whilst this may be true the presence 
of these distributions may actually be attributed to the 
manner in which the networks undergo construction. 

The mechanism of preferential attachment Q has been 
used successfully to model the development of many 
non-biological networks 0, 0], reproducing power laws 
with the correct range of exponents. For biological sys- 
tems, the process of node duplication is more appropri- 
ate and models of this type have again been efficacious 
0, 0, 0, H, 0, 0| . These two mechanisms have a level 
of equivalence though and a common theme amongst 
their application is the reliance upon node growth in the 
dynamics. 

Approaches using fixed node numbers or node dele- 
tions have been proposed that use a variety of mecha- 
nisms to achieve scale-free and other distributions 0, 
I131 113 . nil llij . In this letter we propose a simple 
model of duplication and deletion where the node num- 
ber remains constant throughout the network evolution. 
The duplication process can be incurred perfectly or with 
a relaxation in the fidelity in order to represent partially 
correlated duplications. 

With perfect duphcation our model produces a 1/fc 
power law degree distribution in the limit of vanishing 
connectance. An ansatz for the functional form of the 
degree distribution is proposed based upon symmetry 
considerations and is fitted well by the simulation data. 
With certain approximations we provide a solution to 
the mean field degree evolution equation which supports 
the presence of an exponent, 7=1. When we relax the 
condition of perfect duplication we see distributions bet- 
ter described by the exponential form. Ecological field 
data often displays such distributions [itI and a pre- 
vious model of species interaction networks has presented 
them also ^ . It could be argued that ecological systems 
develop via processes of extinction and correlated specia- 



tion with steady state species numbers so the appearance 
of similar degree distributions here is of relevance. 

The model - We consider a network of N nodes where 
edges are undirected and self-connection is excluded. The 
initialised edge set is transient under the proposed dy- 
namics so can be taken from any arbitrary distribution. 

The edge evolution progresses via non-local updates 
but is nevertheless very simple. At a timestep, we ran- 
domly select with uniform probability a node and all 
its associated edges to be deleted from the network. A 
second parent node is then randomly selected from the 
remaining A^— 1 nodes and is duplicated in the form 
of a daughter node. All nodes connected to the par- 
ent are now given connections to the daughter with 
probability Pe- AH nodes not connected to the par- 
ent are given connections to the daughter with proba- 
bility, P„. The determination of an edge between the 
daughter and the parent is made with probability. Pp. 
As will be shown, at mean field level the system ad- 
mits a single attracting fixed point in the connectance, 
C— >Co, which is determined by the duplication probabil- 
ities (EqO. By using this relation we select the prob- 
abihties, Pe,Pn, dependently such that Pp—Co, so al- 
lowing FeG[Co, 1] to be used as a control parameter for 




FIG. 1: Degree distributions for iV = 200, Co = 0.01 systems 
produced using the imperfect duplication process. From short 
to long tail we have Pe={0.01, 0.25, 0.75, 0.95, 0.99, 0.999}. 
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FIG. 2: Perfect duplication degree distributions for networks 
with N={25, 50, 100, 200, 400, 800, 1600}, Co=0.01. Each has 
been normalised with the exclusion of the fe=0 support and 
rescaled to overlap for visual purposes. 



the fidelity. As a result, the dynamics range from fuUy 
random, Pf,=Pn=Pp=Co, generating binomial networks, 
to fully correlated, Pe—l,Pn=0,Pp—Co, where daughter 
and parent are indistinguishable. 

The ensemble averaged degree distributions for various 
fidelities are shown in Fig.(^. For the uncorrelated, ran- 
dom duplication process, Pe=P„=Co we see binomially 
distributed networks as expected. As we increase the fi- 
delity, the distributions exhibit longer tails and conform 
more to exponential than they do binomial. This pro- 
gression continues until we have a process of perfect du- 
pHcation at which point the distributions become power- 
law-like. 

Shown in Fig. (01, are the ensemble averaged degree 
distributions for perfect dupUcation networks with Cq = 
0.01 over increasing system sizes. For finite but small 
values of the connectance the distributions achieve forms 
that conform closely to a power law. As we reduce Co to- 
wards zero, the functional forms become strongly power- 
law-like with exponents, 7~1. The networks achieved 
here are extremely sparse and the majority of the sup- 
port is held in the P(fc=0) degree. So the form we see 
here represents (near) scale-free fluctuations above the 
k=0 vacuum. In the limit that Cq—^O a finite sized sys- 
tem asymptotically achieves the absorbing state of zero 
connectance. As we take the thermodynamic limit in the 
network size, though, we have a system where P{k—0) 
tends to unity and the remaining distribution achieves 
power-law form but with the support for P(fc>0) tend- 
ing to zero. The simulations suggest that in these limits 
we have a network of zero connectance with 1 /k scale- free 
fluctuations in the degree. It is also interesting to observe 
a mirror symmetry between distributions of complemen- 
tary connectance {Cq'—I—Cq), in Fig.©. The distribu- 
tion for Co=0.49 is also shown, demonstrating their grad- 
ual convergence to the same uniform profile at Co=0.5. 

Symmetries - We propose a functional form for the 



degree distributions acquired with the perfect dupUca- 
tion process, based upon symmetry requirements from 
both theoretical and empirical information. Note that 
the probability of duplicating an edge, Pe—l, whilst the 
probability of duphcating a lack of edge is, 1— P„ = l. 
Additionally, the parent-daughter edge is created with 
probability, Pp, and complementary to this a non-edge 
is, 1— Pp. There are no other aspects to consider in this 
process. Therefore, if we shift our frame of reference 
from edge-centric to non-edge-centric we may state that 
the degree distribution has the following symmetry, 

Pik';Pp,N) = P{l-k';l~Pp,N), (1) 

where we have used a scaled variable for the degree, 
k'=j^^. This duality relation is quite general with mul- 
tiple classes of solutions but we can constrain the pos- 
sibilities by using the empirically observed mirror sym- 
metry. When comparing distributions of complementary 
connectance Cq'—I—Cq (Fig|2J), they appear to be mirror 
symmetric versions of one another on the log-log scale, a 
feature that can be represented by the relation, 

logP{k'; Pp, N) = -logPik': 1-Pp, 7V)+A(Pp, N), (2) 

where A{Pp,N) is a constant for a given Pp and N. If 
we use the previous symmetry of Eq.l^, this may be 
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FIG. 3: Perfect duplication degree distributions for A^=200, 
Co=0.01,0.99 (top); Co=0.49 (bottom). Curve fits are ap- 
plied using Eq.(0J. 
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reformulated as, 

P(fc'; Fp, N)P{l~k'- Pp, N) = A{Pp, N). (3) 

So with these two symmetries in mind we now propose 
an ansatz for the degree distribution of the form, 



'P(0;Pp,7V), 
Pik';Pp,N) = \G{Pp,N)[- 
P{l;Pp,N), 



k' 



-{l-2Pp 



fc'=0; 

0<fc'<l; 

k'=l. 



(4) 

This is consistent with both the theoretical and observed 
symmetries and agrees strongly with the distributions 
acquired from the simulations (FiglS))- For the uniform 
degree distribution of Co=0-5 our functional form also 
fits well. The exponent of the function tends to zero at 
this point leaving P{k'; 0.5, N)=constant as required. 

Mean field theory - To demonstrate the existence of an 
attracting fixed point in the connectance we now consider 
the mean field dynamics of the number of edges, Et , 



= Et-kt+Pekt+l+Pn{N-2~kt+l)+Pp, 

2 



t+2 — 

kt — , 

N ' 



N-1 



{Et-kt) 



(5) 



The second term accounts for the loss of edges when a 
node with mean degree, kt is deleted in stage one. The 
third and fourth terms then represent the increase in 
edges incurred at stage two when one of the remaining 
iV— 1 nodes is duplicated. Imperfect duplication is made 
possible at this point by the inclusion of the probabilities, 
Pe and P„. The fifth term represents the contribution 
from the daughter-parent connection. 

As we only observe the system after every two steps we 
may rescale the time variable increment, (i+t')— 
Thus Eq.© may be expressed as a simple recursion re- 
lation which admits a single attracting fixed point that 
can be expressed in terms of the network connectance, 
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P,{N-2)+Pj, 



N-l-{P,-Pn){N-2)- 



(6) 



We now consider the evolution equation for the aver- 
age number of nodes of degree k, Uk, constrained by the 
condition, nk=N, 

nk{t+l)=nk{t)+rR{N,k,t)+TDuiN~l,k,t). (7) 

Here the probabilities, rji{N,k,t) and ri:)u{N—l,k,t) 
concern the effects of removal (R) and duplication (Du) 
respectively. They consist of the following contributions 
(we suppress for convenience the timestep variable t and 
the size of the network, TV or N—1): 



Tnik) 



T},{k)+T%ik+l)~T%{k) 
^uik) 



+^Duik~l) 



^ Dii 



ik) 



(8) 



(9) 



The separate terms represent the effects of the pro- 
cess on the distinct node types, indexed as: removed (r), 
adjacent (a), parent (p), daughter (d). 

The direct effect on rik of removing a node of degree, 
k is to decrease its support by one. The probability of 
selecting a node of degree A: is and therefore, 



(10) 



Upon removal, the degree of the nodes connected to the 
removed node will decrease by one. For this we need 
the Edge probability, PEd{ki,k2,q), the probability that 
a node of degree ki is connected to q nodes of degree fc2. 
This will be estimated below fEa ll7(l . Here we write, 



N 
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q=l 
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The first sum is over the degree of the removed node, 
the second sum is over the number of nodes of degree 
fc = 0, 1, A^— 1 the removed node is connected to. 

When a daughter is introduced it receives an edge to 
the parent with probability Pp. Thus the parent increases 
its degree by one. A node of degree k is selected for 
duplication with probability nk/{N—l), 



N-1' 



(12) 



The daughter can add support to Uk by an amount given 
by the probability of finishing with k edges, 



TiJk) = PpK{k~l) + {l-Pp)K{k). 



(13) 



where, 



K{k) = Probjdaughter receives k edges to nodes 
different from the parent} 

N-2 m[fcp,fe] m[N ~2-kp,k\ 



(14) 



EE E 

fep=0 fci=0 



fe2=0 



A-1 



5{ki-\-k2—k)^{ki,k2, kp 



where m in the limits represents, minimum. This expres- 
sion corresponds to when the daughter inherits fci edges 
associated with nodes connected to the parent with each 
attachment occurring with probability, Pg. The remain- 
ing k2=k—ki edges are established with probability, P„ 
and are associated with nodes not connected to the par- 
ent. The factor ^{ki,k2,kp) is the probability that ki 
edges derive from the kp nodes connected to the parent 
and k2 from those not connected to the parent. 



n{k_ 



N~2~k 
k2 



V \ P^^(^l — P 'jN-2-kp-k2 



(15) 
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The duplication process also effects potential adjacent 
nodes connected or unconnected to the parent. A par- 
ent of degree kp is connected to q nodes of degree k with 
probability Psdikp, k, q). These q nodes will as a result of 
the duplication with probability Pg become degree fc+l 
nodes. There are N~2~kp nodes the parent has no link 
to. Among these nodes q will receive a new edge to the 
daughter and become degree fc+l nodes with probabil- 
ity PnPEd{N~2'-kp, k, q). The total contribution to the 
transition probability is, 

^nu(k) = E E E E EE ^('^1+^-'^) 

g— kp—OK,i—0 K.2— 91—0^2—0 

q^PEd{kp,k,K,) 2 ) P^Hl-Per-"' 

PEd{N-2-kp,k,K2) 2 ) PTi'^-Pur'^"'- 

(16) 

We estimate PEdiki,k2,q) by using the following urn- 
model consideration. We have as many particles in the 
urn as there are edges, i.e, M= J2k knk- We imagine that 
the particles are of two types, A and B. The former cor- 
responds to \A\—k2{nk2—Ski,k2)'^ the edges connecting to 
nodes of degree ^2. The latter to |i3|=M— 1^|; the edges 
leading to nodes with degree different from ^2. We pick 
ki particles from the urn. These are the edges which will 
be connected to the single node of degree ki under consid- 
eration. The probability that among these ki particles, 
q are of type A and ki—q are of type B is, 



PEd{ki,k2,q) = 



\A\_. 
M' 



l- 



AP 



(17) 



It is straightforward to check that these equations sup- 
port the normalised stationary solution P{k) oc 1/fc in the 
double limit Pn—^0 and P^— >1, when only leading contri- 
butions in the sum over Psd from Ea. H16|) are included. 
This implies the relevant correlations make PEd{ki, k2, q) 
decay rapidly to zero with increasing ki and q. Measured 
degree-degree correlations in the numerics support this. 

Conclusion - We have shown that a power-law degree 
distribution may be achieved without the usual require- 
ment of node growth. The ansatz for the degree distri- 
bution suggests that the associated exponent is equal to 
unity, and the solution to the degree evolution equation 
in the required limits supports this. Both the processes of 
preferential attachment and duplication can give a range 
of exponents but achieving an exponent this low is un- 
usual. Biological systems and in particular ecological sys- 
tems often exhibit low range exponents when power laws 
occur, 7 ~ (1,2) [ll, so it is interesting to see this re- 
produced in our model. Ecological systems often present 
degree distributional forms other than the power-law and 
our model shows that partially correlated duplication can 



lead to distributions better described as exponential. The 
appearance of non-binomial distributions in ecology may 
just be an inherent consequence of the system dynamics. 

The duplication process in growth models naturally in- 
corporates preferential attachment which allows a strong 
comparison between the two approaches. But for our 
model the equivalence between the two mechanisms be- 
comes blurred. By utilising deletion as well as duplica- 
tion, higher degree nodes again have greater probabili- 
ties of increasing their degree during duplication events. 
But conversely they have greater probabilities of reduc- 
ing their degree during the deletion events ie. preferential 
detachment. The 'more begets more' adage seems inap- 
propriate here as more also begets less. 
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